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NEUTRON STARS WITH A STABLE, LIGHT SUPERS YMMETRIC BARYON 
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ABSTRACT 

If a light gluino exists, the lightest gluino-containing baryon, the S° , is a possible candidate for self-interacting 
dark matter. In this scenario, the simplest explanation for the observed ratio Q<jm/^& ~ 6 — 10 is that 77150 ~ 
900 MeVc~ 2 ; this is not at present excluded by particle physics. Such an S° could be present in neutron stars, 
with hyperon formation serving as an intermediate stage. We calculate equilibrium compositions and equation of 
state for high density matter with the S° , and find that for a wide range of parameters the properties of neutron 
stars with the S° are consistent with observations. In particular, the maximum mass of a nonrotating star is 
1.7—1.8 M Q , and the presence of the S° is helpful in reconciling observed cooling rates with hyperon formation. 

Subject headings: dark matter-dense matter-elementary particles-equation of state-stars: neutron 



1. INTRODUCTION 

The interiors of neutron stars offer a unique meeting point 
between astrophysics on one hand and nuclear and particle 
physics on the other. The macroscopic properties of neutron 
stars, such as mass, radius, spin rate and thermal history de- 
pend on the microscopic nature of matter at very high densi- 
ties. Neutrons stars are, in this sense, cosmological laboratories 
of physics at densities exceeding the nuclear saturation density 
(po ~ 0.16 fm~ 3 ). One aspect of matter at supernuclear densi- 
ties is the possible appearance of hadronic degrees of freedom 
beyond nucleons. The principle candidates are hyperons, me- 
son condensates and a deconfined quark phase (see, e.g., Glen- 
denning 1996, Prakash et al. 1997 for reviews). 

In this work we wish to illuminate another possibility, the 
appearance of an i?-baryon. We are motivated by the obser- 
vation of Spergel & Steinhardt (2000) that the observed prop- 
erties of dark matter (DM) on scales of order 10-100 kpc are 
better described by a theory with self-interacting DM than by 
conventional, weakly interacting DM. Since the required cross- 
section to mass ratio, crxx I mx ~ lb/GeV, is characteristic of 
hadronic physics at low energies, a light stable neutral hadron 
would be a welcome feature of physics-beyond-the standard 
model, and such particles arise naturally in theories with a very 
light gluino. In some sense, this is the minimal extension of the 
standard model needed to produce the required DM particle, 
and therefore merits examination. 

A supersymmetric version of the standard model contains 
a massless gluino (g), the spin- 1/2 superpartner of the gluon; 
therefore it is endowed with a corresponding chiral symmetry 
called i?-invariance. When supersymmetry is broken, the R- 
invariance may or may not break. If it does not, the gluino 
mass arises only from loops sensitive to the Higgs vevs (which 
do break i?-invariance) and is naturally of order 100 MeVc~ 2 . 
Although considerable effort has been devoted to excluding it, 
a very light gluino remains a possibility - mainly because of 
the difficult of making firm, first-principles predictions for non- 
perturbative phenomena in QCD. The question deserves very 
serious scrutiny because the existence of a light gluino would 
resolve several puzzles. These issues and the constraints im- 
posed by experiment are reviewed by Farrar (1998, 2000). See 
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also Farrar & Gabadadze (2000) and references therein. 

The lightest supersymmetric "i?-hadrons" are the spin- 1/2 
gg bound state (R° or glueballino) and the spin-0 baryon 
udsg bound state, the S° . The S° is thus a boson with baryon 
number unity and a strangeness charge — 1. The strong color 
hyperfine interaction unique to the udsg in a flavor singlet state 
makes the mass of the S° much less than that of all other R- 
baryons (Farrar 1984; Buccella & Farrar 1985), and probably 
less than 2 GeVc -2 , which would make it stable through con- 
servation of R— parity (Farrar 1996). 

In the very early universe the S° remains in equilibrium with 
nucleons via processes such as R a N <-> S !^, down to T rj 35 
MeV (Farrar, Spergel & Steinhardt 2000; see also in Wandelt 
et al. 2000). Hence for a suitable S° mass, slightly less than the 
proton mass, standard freezeout estimates lead to a density ratio 
ft s o /f2s « 6 — 10, as required if the S° is to be the dark matter. 
To posit a new baryon lighter than a nucleon is quite shocking! 
But it is known that the color hyperfine force is very strong in 
this channel, and there are innumerable examples of the inad- 
equacy of phenomenological models of non-perturbative QCD 
phenomena. The possibility of such a light S° should not be dis- 
missed on theoretical grounds, and we must seek observational 
constraints. One possible astrophysical constraint is whether 
this new baryon can form in neutron stars, and if so, would it 
imply unacceptable neutron star properties? 

Nucleons cannot convert directly into the S° , even if the 
latter is less massive. Conservation of R— parity requires that 
i?-hadrons are produced in pairs, forbidding direct conversion 
of individual nucleons into S° particles. The weak neutral cur- 
rent process nn — » S°S° requires simultaneous creation of two 
units of strangeness (AS — —2). Strangeness-changing neu- 
tral currents (AS = ±1, AQ = 0) are highly suppressed in the 
standard model due to the GIM mechanism. Thus, nn — » S°S° 
is doubly GIM suppressed; crude estimates indicate that its 
rate is of order the inverse Hubble time (Farrar et al. 2000). 
However processes such as A A — > S°S° (AS* = 0) and even 
T,~p — > S°S° (charged current, AS = -1) are not GIM sup- 
pressed, so that if the conditions in interiors of neutron stars 
allow for hyperons to appear (as is predicted by many works 
on this issue, see, e.g., Glendenning 1985, Weber & Weigel 
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1989, Knorren, Prakash & Ellis 1995, Schaffner & Mishustin 
1996, Balberg & Gal 1997, Huber et al. 1998, Baldo, Burgio 
& Schulze 2000), these hyperons could serve as a doorway for 
S° formation, which might affect observable properties of neu- 
tron stars. 

Our goal in this letter is to explore whether a light S° can 
form in neutron stars and assess whether the resulting global 
properties of these stars are compatible with observations. In 
§ 2 we extend an equation of state for high density matter with 
hyperons to include the S° . The equilibrium composition with 
the S° is found in § 3, followed by a discussion of the corre- 
sponding properties of neutron stars in § 4. 

2. S° INTERACTIONS IN HIGH DENSITY MATTER 

In order to include S° particles in an equation of state of high 
density matter, we must specify their interactions with other 
baryons and among themselves. Here we do so in the frame- 
work of an effective equation of state (Balberg & Gal 1997, 
Balberg, Lichtenstadt & Cook 1999), where the potential felt 
by a single, isospin singlet particle of species x placed in mat- 
ter composed of species y at a baryon number density p y is 



Vx(Py) = a X yPy + C X yPl + W Py 



(1) 



In equation (1) a is negative and represents a long range, attrac- 
tive component in the interaction, c is positive, corresponding to 
an intermediate range repulsive component, and w is also pos- 
itive, representing a universal repulsive short range interaction, 
presumably arising from the quark structure of the baryons. The 
powers 7 and 9 must be greater than unity, so that repulsion 
dominates at high densities. 

For nucleons and hyperons, the values of the coefficients can 
be constrained by experimental data. Since no such data exists 
for the S° , we resort to phenomenological arguments to esti- 
mate the magnitude of these terms. The attractive component 
is generally associated with two-pion or two-kaon exchanges. 
This interaction is strong because the intermediate state in the 
s-channel can be the same as the initial state, hence the ampli- 
tude has a small energy denominator. However the flavor sin- 
glet S° cannot exchange tt or K without becoming a much more 
massive flavor octet i?-baryon (Buccella & Farrar 1985). Thus 
the long-range attractive interaction is strongly suppressed for 
the S° . A more empirical way to arrive at the same conclusion 
is to note that if there were a significant long-range attractive 
force the S° would bind to nuclei (and to itself), and nuclei 
with exotic mass values have not been observed. We there- 
fore assume that the attractive component is negligible for the 
S° —related interactions. The intermediate repulsive component 
may be associated with the exchange of vector mesons, such as 
the p and r\. Only the flavor-singlet component of these vector- 
mesons can couple to the S° , so we roughly assume that the 
magnitude of this term (the value of c in Eq. [1]) is half of its 
value for the A. The universal short-range interaction for the 
S° is taken as identical to that of the baryons. Here we use EoS 



1 of Balberg et al. (1999), where 7 = \ and 6 



7, so that for 



the S° a = and c = 329.5 MeVfm 4 for all its interactions, 
and w = 223.6 MeVfm 5 . Below we refer to this choice of pa- 
rameters along with an S° mass of rag = 900 MeVc~ 2 as our 
nominal set. Most of our conclusions prove insensitive to these 
choices, as will be seen below. 



3. COMPOSITION OF NEUTRON STAR MATTER WITH THE S° 

The composition of neutron star matter is determined 
by chemical equilibrium for all weak-interaction processes. 
Specifically, these are B\ — > Bitvg, where B\ and Bi are 
baryons, I is a lepton and ut is its corresponding neutrino (or 
antineutrino). Neutron star matter is well approximated as cold 
(T = 0) and neutrino free (fj, v = 0) so that all equilibrium 
conditions are summarized concisely in the form: 



(2) 



where pB t is the chemical potential of the baryon species Bi, 
fi n and p e are the neutron and electron chemical potentials, re- 
spectively, and is the electric charge of the baryon Bi. The 
equilibrium composition at any baryon number density, ps, is 
then determined by two additional conditions: T,iXBi = 1 and 
charge neutrality, HiXBSBt + ^i x eqe = 0, where and 
xi denote the fraction per baryon for each baryon and lepton 
species, respectively. 

Equation (2) does not fully apply to the case of the S° , since 
the process nn — ► S°S° is suppressed, even if energetically fa- 
vorable. The S° must maintain a direct equilibrium only with 
hyperons, i.e., 2/igo = 2/iyo, and 2/igo = p Y - + AV but 
2/igo = 2/i„ is not required (where Y a and Y~ denote a neutral 
and a negatively charged hyperon, respectively). The signifi- 
cance of this distinction is that a combination of S° and nucle- 
ons may coexist even while p s o < \i n , as long as no hyperons 
are present. On the other hand, a finite hyperon fraction will 
indirectly impose that /igo = ji n , so then equation (2) applies 
to all species. 

We emphasize that even if /igo < p n , nuclear matter can- 
not "burn" exothermally into S° matter. Unlike the case of 
"strange-quark-stars" (Witten 1984; Alcock, Farhi & Olinto 

1986) , where a deconfining front converts an entire neutron star 
into the (hypothetically) lower energy strange quark matter if 
a seed of of strange quark matter is somehow formed (Olinto 

1987) , there exists no equivalent stimulant to initiate conver- 
sion of nuclear-matter into S° -matter, regardless of density or 
composition. 

In Fig. 1 a we show for reference the equilibrium composition 
of neutron star matter as a function of baryon density, using 
the "standard" equilibrium composition of matter with nucle- 
ons and hyperons but no S° (hereafter, NH matter) based on 
EoS 1 of Balberg et al. (1999), assuming that the equation of 
state is valid up to densities as large as 10/?o- The results, with 
hyperon formation beginning at about 2/?o, are typical of mat- 
ter with hyperons (Glendenning 1985; Weber & Weigel 1989; 
Knorren et al. 1995; Schaffner & Mishustin 1996; Balberg & 
Gal 1997; Huber at al. 1998; Baldo et al. 2000). 

Figure lb depicts the equilibrium composition when the S° is 
included according to the assumptions of § 2. Due to the 
low mass of the S° , its chemical potential is lower than that 
of all hyperons, even when the fractions of the latter are in- 
finitesimally small. As a result, all hyperons that form convert 
to S° particles, and the matter consists only of nucleons and 
S° . The composition is determined by the nucleons being in 
equilibrium with an exactly zero fraction of hyperons 4 . Even 
though absent from the matter, the hyperons do play a role in 
determining the amount of S° particles which are created. 



4 In reality, some additional nucleons would convert to hyperons through quantum and thermal (T ~ O.OlMeV in neutron star cores) fluctuations, 
and these hyperons would decay to S° 's. The actual equilibrium composition would be one where the chemical potential of the nucleons is slightly 
below that of zero-fraction hyperons at every density (in other words, fluctuations would deplete a thin outer layer in phase space of nucleons). 
Assuming that the probability of fluctuations decreases exponentially with the energy gap, the quantitative effect on the results should be small. 
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FIG. 1. — Equilibrium compositions for high density matter for EoS 1 
of ref. [1]: (a) matter with hyperons when the S° is ignored, (b) includ- 
ing the S° with m s o = 900 MeVc~ 2 , and (c) same as (b) but with 
m s o = 975 MeVc- 2 . 

Since the S° is never in chemical equilibrium with the nu- 
cleons (or the zero fraction hyperons), the equilibrium compo- 
sition would be identical for any ms < 900 MeVc~ 2 . In 
fact, for the nominal interactions, the same composition is also 
found for mg < m p . A qualitative change is found when 
ms is larger: in Fig. lc we show the equilibrium composition 
for nominal S° interactions, but with m s o = 975 MeVc~ 2 . 
At low and intermediate densities the composition is identical 
to that of the mgo = 900 MeVc~ 2 case, but beginning at a 
density of 0.9 fm -3 the larger mass of the S° makes it favor- 
able for a finite fraction of hyperons to coexist in the matter. 
The general equilibrium condition (Eq. [2]) is now satisfied for 
all particle species. Owing to the S° having only repulsive in- 
teractions, the hyperons achieve chemical equilibrium with it 
even though it has a lower mass and no kinetic energy. More- 
over, for larger densities, it becomes increasingly favorable to 
reduce the S° fraction as the equilibrium composition gradu- 
ally approaches that of NH matter. With similar calculations we 
find that for nominal interactions, if m s o ~ 1030 MeVc~ 2 the 
S° appears only after some hyperons already exist in the matter, 
and if m s a > 1060 MeVc~ 2 S° formation is suppressed for the 
entire range of densities. 

4. EQUATION OF STATE AND MASS SEQUENCES 

We now turn to examine the equation of state (pressure 
vs. mass-energy density, hereafter EoS) of high density matter 
with the S° . These are shown in Fig. 2, including the models 
used for Fig. 1 and also a model labeled N for nuclear matter, 
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FIG. 2. — Equation of state for matter of the compositions shown in Fig. 1. 
Matter with the S° : nominal values (see text) - thick solid line; nominal inter- 
actions and m s o = 975 MeVc~ 2 - dashed line; nominal mass and increased 
repulsion - dot-dashed line. Also shown are the equations of state for nucleon 
and hyperon matter - N+H, and nuclear matter - N. 

where hyperon (and therefore also S° ) formation is ignored. 

The hierarchy of EoSs in Fig. 2 is mostly determined by the 
number of available species in each model. Nuclear matter nat- 
urally has the stiffest EoS, followed by matter with the S° . 
Even though the S° , being a boson, does not exert any Fermi 
pressure, the fact that four species of hyperons are suppressed 
due to its appearance makes such matter stiffer than NH mat- 
ter, which is clearly the softest. An exception occurs just above 
the threshold density for hyperon formation, where NH mat- 
ter includes only a single hyperon species (the so that the 
contribution of its Fermi pressure does make NH matter stiffer. 

Although the mass of the S° does not affect its contribution 
to the pressure, it does contribute to the mass energy-density. 
Hence a lower mso yields a stiffer EoS, while a larger value 
yields a softer one. For m s o = 975 MeVc -2 such softening 
is enhanced at a density of ~ 1.8 x 10 15 gmcm~ 3 , when hy- 
perons begin to accumulate; as the S° fraction is depleted this 
EoS gradually approaches that of NH matter. It is noteworthy 
that increasing the repulsive power of the S° does not lead to 
significant stiffening of the EoS. We demonstrate this with an- 
other model shown in Fig. 2 when m s o = 900 MeVc~ 2 but 
c = 494.5 MeV fm 4 (50% larger than the nominal value). 
Magnified repulsion leads to an increase of the S° chemical po- 
tential, which in turn enhances hyperon formation. The result- 
ing composition of this model is quite similar to that of Fig. lc, 
and the EoS, while slightly stiffer than nominal until hyperons 
appear, is then softened. We emphasize that this is a mani- 
festation of the ability of hyperons to exert "pressure control" 
(Balberg & Gal 1997; Balberg et al. 1999), since the extent of 
their formation tracks the dependence of particle interactions 
on density. 

Given an EoS, the structure of a spherical, non-rotating 
neutron star is found by numerically integrating the Tolman- 
Oppenheimer-Volkoff equations of general relativistic hydro- 
static equilibrium (see, e.g., Misner, Thorne & Wheeler 1973). 
The resulting mass sequences (stellar mass as a function of the 
assumed central density) for the EoSs of Fig. 2 are given in 
Fig. 3. All sequences are identical up to about 1.1 M Q , where 
the mass sequences diverge according to the hierarchy of EoSs. 
The maximum mass found for matter with the S° is 1.81 M 
for the nominal case, 1.68 M Q for m s o = 975 MeVc~ 2 , 
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FIG. 3. — Nonrotating neutron star sequences for the five equations of state 
shown in Fig. 2. 

and 1.78 M Q for m s o = 900 MeVc" 2 but enhanced S° re- 
pulsive interactions. In both latter cases, hyperons coexist with 
the S° only in the largest mass stars, those of mass closer than 
0.1 M Q to the maximum mass in each case. We note that for 
EoS 1, the maximum masses for nuclear matter and NH matter 
are ^ 2.18 and ~ 1.57 M Q respectively. 

At present, the primary observational constraint on high den- 
sity equations of state is that the maximum mass must be at 
least 1.44 M , the inferred mass of PSR 1913+16. All EoSs 
discussed above, with and without the S° , are acceptable. This 
limit does allow to rule out very soft EoSs, such as would arise 
if the S° had a low mass and an attractive component in its in- 
teractions. Such interactions would make a pair of S° similar to 
the hypothetical H dibaryon - a six quark flavor singlet, boson 
(Jaffe 1977a, 1977b) - whose interactions are supposed to be 
similar to those of two A hyperons combined. The EoS for mat- 
ter with the H would indeed be too soft to sustain a 1.44 M Q 
star, unless its mass is significantly above 2m\ (Glendenning 
& Schaffner-Bielich 1998). The expected strongly repulsive in- 
teractions of the S° are thus crucial for allowing its existence in 
neutron stars. 

5. CONCLUSIONS AND DISCUSSION 

Our primary conclusion is that the existence of a light 
( m S° ^5 m p) S° particle is consistent with the fundamen- 



tal observed properties of neutron stars, namely, their inferred 
masses. This conclusion applies for a wide range of models 
for high density matter with the S° . In particular, the maxi- 
mum mass of a static neutron star under nominal assumptions 
regarding the S° is found to be about 1.8 M Q . We note that 
since it suppresses hyperon accumulation, the appearance of the 
S° yields a stiffer equation of state and higher maximum mass 
than obtained with nucleons and hyperons without S° . There- 
fore, identification of a relatively high mass neutron star would 
be easier to reconcile with an S° component than in its absence 
(where hyperons accumulate undisturbed). Some indication for 
the existence of such massive neutron stars is implied by re- 
cent analyses, although uncertainties are still significant (van 
Kerkwijk 2000). We emphasize that current data on neutron 
star masses is not restrictive enough to decisively constrain the 
properties of high density matter. 

The presence of S° particles could have a significant and wel- 
come effect on the cooling rate of neutron stars, if they do in- 
deed replace hyperons in the cores. Hyperons can participate in 
the direct Urea cooling process (Bi — > B 2 ev) while the S° can- 
not. The threshold fraction required for hyperons to drive direct 
Urea cooling processes is of the order of 0.01 (Prakash et al. 
1992), so neutron stars with hyperons are susceptible to rapid 
cooling, which seems inconsistent with observations (Schaab 
et al. 1996). Hyperon direct Urea cooling can be partially sup- 
pressed if the A hyperons pair to a 1 So superfluid state (Schaab, 
Balberg & Schaffner-Bielich 1998), but the extent of suppres- 
sion is sensitive to the details of the equilibrium composition 
and gap energy (Page et al. 2000). On the other hand, S° forma- 
tion would prevent the hyperon direct Urea, or restrict it only 
to the highest mass stars, so that very moderate cooling rates 
are attainable. Future measurements of surface thermal emis- 
sion from young neutron stars with the CHANDRA and XMM 
satellites may lead to stronger constraints on cooling rates; if so 
this will yield valuable insight into the internal structure of high 
density matter. 

The research of GRF is supported in part by NSF-PHY-99- 
96173. The idea that the S° could be a dark matter candidate, 
and constraints on such a scenario, were developed by GRF in 
collaboration with D. Spergel and P. Steinhardt; that work will 
be reported in greater depth elsewhere (Farrar et al. 2000); a 
brief description can be found in (Wandelt et al. 2000). SB and 
TP thank NYU for hospitality which facilitated this research. 
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